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Intragenic integration in DLC1 sustains factor VIII expression in 
primary human cells without insertional oncogenicity 

J Sivalingam 1 ' 2 , TT Phan 3 ' 4 and OL Kon 1 ' 2 

Techniques enabling precise genome modifications enhance the safety of gene-based therapy. DLC1 is a hot spot for phiC31 
integrase-mediated transgene integration in vitro and in vivo. Here we show that integration of a coagulation factor VIII transgene 
into intron 7 of DLC1 supports durable expression of factor VIII in primary human umbilical cord-lining epithelial cells. Oligoclonal 
cells with factor VIII transgene integrated in DLC1 did not have altered expression of DLC1 or neighbouring genes within a 1-Mb 
interval. Only 1.9% of all expressed genes were transcriptionally altered; most were downregulated and mapped to cell cycle and 
DNA repair pathways. DZ_C7-integrated cells were not tumourigenic in vivo and were normal by high-resolution genomic DNA copy 
number analysis. Our data identify DLC1 as a locus for durable transgene expression that does not incur features of insertional 
oncogenesis, thus expanding options for developing ex vivo cell therapy mediated by site-specific integration methods. 
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INTRODUCTION 

One of the fundamental requirements for developing gene-based 
cell therapy for disorders caused by deficient production of 
specific proteins, such as the haemophilias, is durable expression 
of the corrective transgene product. This can be optimally 
achieved by stably integrating the appropriate transgene into 
the genome. Several clinical trials that used integrating viral 
vectors have provided proof-of-principle by successfully correcting 
the disease phenotype of several monogenic disorders. 1-6 
However, serious adverse outcomes that later emerged in some 
of these studies brought to light oncogenic complications 
incurred by random integration of transgenes. It is now 
understood that transactivation of neighbouring oncogenes 
such as LM02, CCND2, BMI1 and the MDS-EVI1 complex 
following retroviral vector integration can lead to clonal cell 
expansion, myelodysplasia and overt leukaemia. 7-9 These 
oncogenic complications motivate the continuing search for 
non-viral methods that integrate therapeutic transgenes in safe 
genomic regions. 10 

The non-viral phiC31 integrase system induces stable expres- 
sion of transgenes that are integrated into endogenous pseudo 
attP sites in mammalian genomes. 11 " 14 It has been successfully 
used to correct deficiencies of fumarylacetoacetate hydrolase, 15 
factor IX, 16 factor VIII 17,18 and dystrophin 19 in murine disease 
models. Bioinformatic analysis predicted that phiC31 integrase 
could potentially mediate integrations into about 370 different 
genomic sites. 11 In practice, however, integrations have been 
found experimentally in only a small subset of these sites. 
Along with others, we have identified a few sites in the 
human genome where transgenes integrate with high 
frequency. Among these hot spots are hitherto uncharacterized 
loci in 8p22 (refs 12 < 13 < 17 < 20 > and 19q 13.31. 11 

We have reported that phiC31 integrase-modified primary 
human cord-lining epithelial cells (CLECs) durably expressed a 



human coagulation factor FVIII (FVIII) transgene and corrected the 
disease phenotype when implanted in FVIII-deficient mice. 17 Up to 
40% of FVIII transgene integrations in a mixed population of CLECs 
occurred in 8p22, and clonal cultures of genome-modified CLECs 
showed no clear markers of genotoxic risk. This drew our attention 
to the 8p22 locus that appeared to be both permissive for durable 
transgene expression and potentially safe against oncogenic risk. 
Here, we show that transgene integration into intron 7 of DLC1 in 
8p22 leaves a minimal and benign footprint in the genome and 
transcriptome, and does not induce tumourigenic behaviour in 
genome-modified cells. Our data identify this site as a likely safe 
harbour for gene-based cell therapies that require integration of 
transgenes. 

RESULTS 

A high expression human-porcine FVIII transgene 
We designed and assembled a hybrid human-porcine FVIII 
complementary DNA (cDNA) that was similar to a construct 
reported to induce 10-fold higher expression than human FVIII 
cDNA. 21 A plasmid encoding B domain-deleted (BDD) hybrid 
FVIII cDNA comprising porcine A1 and A3 domains, a residual 
human B domain retaining eight glycosylation sites, and 
human A2, C1 and C2 domains. To compare the efficacy of the 
human-porcine hybrid BDD FVIII cDNA with BDD human FVIII 
cDNA F309S, 17 CLECs were co-electroporated with plasmids 
expressing enhanced green fluorescent protein and either 
human or human-porcine FVIII. Both FVIII cDNAs were 
expressed from the human ferritin light chain promoter. After 
normalizing for transfection efficiency, assays of FVIII activity 
secreted by transiently transfected CLECs showed approximately 
5.5-fold higher secretion of human-porcine FVIII (6.13 ± 1.07 mU ml -1 
per 24 h) compared with human FVIII (1.12 ± 0.26 mU ml -1 per 24 h; 
P = 0.003; Figure 1a). 
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PhiC31 integrase-mediated integration of hybrid FVIII cDNA in 
CLECs 

To determine durability of FVIII transgene expression by phiC31 
integrase-mediated genomic integration, CLECs were electropo- 
rated only with paffShybrid FVIII (an c?ff8-bearing plasmid 
encoding the neomycin resistance gene and hybrid FVIII cDNA 
expressed from human light chain ferritin promoter; Figure 1b: 
pctffShybrid FVIII) or co-electroporated with pCMV-Int (a plasmid 
encoding CMV promoter-phiC31 integrase) 22 using the Amaxa 
Nucleofector I device (Lonza, Basel, Switzerland). Electroporated 
CLECs were either unselected (Figure 1b: pctff£hybrid FVIII and 
pctffShybrid FVIII + Integrase) or subjected to G418 selection 
(Figure 1b: paffShybrid FVIII + Integrase G418) to derive cultures 
stably integrated with the pcrtfflhybrid FVIII construct. FVIII 
activity assay (Figure 1b) of conditioned media showed durable 
FVIII expression from G418-selected CLECs that were co-electro- 
porated with pctff£hybrid FVIII and phiC31 integrase (day 6: 
298.30 ± 1 9.01 mU FVIII/1 0 6 cells per 24 h; day 36: 635.6 ± 3.40 mU 
FVIII/10 6 cells per 24 h) but not from CLECs that were similarly 
transfected without G41 8 selection (day 6: 298.30 ± 19.01 mU FVIII/ 
10 6 cells per 24 h; day 36: 8.97 ±2.51 mU FVIII/10 6 cells per 24 h), 
indicating that G418 selection enriched for cells with stable 
genomic integration of the FVIII transgene. CLECs that were 
transfected only with pcrtfShybrid FVIII and not G41 8-selected also 
showed a steep decline in FVIII secretion from day 6 to day 36 (day 
6: 89.02 ± 9.77 mU FVIII/10 6 cells per 24 h; day 36: 7.16±2.40mU 
FVIII/10 6 cells per 24 h), consistent with loss of episomal FVIII 
plasmid with successive cell divisions and negligible genomic 
integration in the absence of the integrase. 

In our earlier study, transgene integrations into 8p22 accounted 
for about 40% of all recovered integration events in a mixed 
population of CLECs. 17 In order to characterize only CLECs with 
8p22 integration at high resolution, we screened stable G418- 
resistant CLECs derived by co-electroporating pcrtfShybrid FVIII 
and phiC31 integrase by junction PCR using primers specific for 
sequences within the FVIII vector and a previously identified 
pseudo attP site in 8p22. 11,20 PCR amplification of both left and 
right integration junctions was evidence for site-specific 
integration of FVIII transgene cassette at 8p22 (Figure 1c). 
Sequencing the integration junction PCR products confirmed 
transgene integration at the 8p22 hot spot. Sequence analysis 
showed a 6-bp deletion in the vector sequence at the left 
integration junction and a 7-bp deletion in the genomic region at 
the right integration junction (Supplementary Figure S1). 



Screening oligoclonal CLECs with 8p22 integration 
Evidence of 8p22 integrations in the bulk-transfected population 
prompted us to derive clonal populations of genome-modified 
CLECs having this specific integration. We examined oligoclonal CLEC 
cultures obtained by flow-sorting four cells into each well of 96-well 
plates for evidence of 8p22-specific transgene integration by direct 
in situ PCR. Using two primer pairs specific to the vector sequence 
and sequence at the 8p22 hot spot, genomic DNA from \n situ lysed 
cells were screened by PCR for the presence of left and right 
integration junctions, indicative of complete transgene integration. 
Intron-spanning primers specific to C1 -domain of FVIII cDNA were 
used to amplify and detect the presence of the integrated transgene. 
A control PCR that amplified an adjacent genomic region in 19q13.42 
served both to verify integrity of genomic DNA and to act as a 
positive control for PCR among all samples. Sixty-six of 72 oligoclonal 
populations were positive for the control PCR amplification. Of these 
positive oligoclones, 16 (24%) amplified both left and right 
integration junctions at 8p22 (Figure 2a). 

Transgene integration in chromosome 8 in oligoclonal CLECs 
was confirmed by fluorescence in situ hybridization (FISH) data, 
which showed close proximity of signals specific to probes for 
chromosome 8 centromere (red signal) and the integrated 
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Figure 1. Enhanced secretion of human-porcine hybrid BDD-FVIII 
and phiC31 integrase-mediated integration by transfected CLECs, 
(a) FVIII activity (Coamatic FVIII assay) in conditioned media of naive 
CLECs or CLECs co-electroporated with an enhanced green 
fluorescent protein reporter gene and either human BDD-FVIII 
F309S cDNA or human-porcine hybrid BDD-FVIII cDNA. Data are 
expressed as mU ml -1 per 24 h normalized to transfection efficiency. 
Each data point shows the mean of triplicate assays of a single 
experiment. An average 5.5-fold increase in FVIII levels (P = 0.003; 
Student's unpaired f-test) was secreted by the hybrid FVIII cDNA 
compared with human FVIII cDNA. (b) FVIII activity of overnight 
conditioned media of three groups of electroporated CLECs: 
(i) pafftfhybrid FVIII, no G418 selection; (ii) paff£hybrid FVIII + 
Integrase, no G418 selection; and (iii) paff£hybrid FVIII + Integrase 
with G418 selection. Conditioned media was collected on day 
6 (before G418 selection) and day 36 (after G418 selection; where 
indicated) post-electroporation. Data are mean±s.e.m.; n = 3 per 
group, (c) Genomic DNA (200 ng) extracted from a bulk population 
of genome-modified CLECs was amplified with primer pairs specific 
for the vector and chromosome 8p22 to detect the presence of left 
and right integration junctions indicative of correct integration at 
the pseudo attP site in 8p22 (done in duplicates), '-ve' denotes 
minus template PCR amplification while ' + ve' denotes amplification 
from genomic DNA isolated from a clonal line of CLEC previously 
identified to have integration at 8p22. Control PCR amplified 
a 900-bp region in 19q13.42 {AAVS1 locus). Amplified products 
were resolved by 1% agarose gel electrophoresis and imaged 
using BioRad Gel Doc 2000 transilluminator and QuantityOne 
software (Bio-Rad Laboratories). 
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transgene (green signal; Figure 2b). This was consistent with 
results of junction PCR (Figure 2a) and amplicon sequencing 
(Supplementary Figures S1c and d). Oligoclonal CLECs with 
transgene integration at 8p22 were further characterized by FISH 
using transgene-specific probes to determine transgene copy 
number (Supplementary Table 1). The majority of cells had a 
single copy of the integrated transgene. However, a small 
proportion of cells (0.6-7%) had two copies suggesting transgene 
integration into other loci besides 8p22. We screened 13 
oligoclonal CLECs derived from the mixed population of G418- 
selected CLECs treated with pcrtfShybrid FVIII + Integrase for FVIII 
expression 40 days post-electroporation. Nine oligoclones 
(Figure 2a; oligoclones 10, 15, 39, 40, 48, 50, 59, 62 and 68) with 
evidence of transgene integration at 8p22 by junction PCR were 
selected, whereas three oligoclones (Figure 2a; oligoclones 18, 30 
and 47) were negative by junction PCR. All nine positive 
oligoclones secreted FVIII at high levels (261 0-5724 mU FVIII/10 6 
cells per 24 h), whereas the three negative oligoclones did not 
secrete FVIII (Figure 2c), showing that transgene integration at 
8p22 sustained durable FVIII secretion by genome-modified cells. 

We proceeded to characterize eight 8p22-modified oligoclones 
for oncogenic and genotoxic alterations, compared with unmo- 
dified CLECs. 

DLC1 expression 

As the transgene integration site in 8p22 is within intron 7 of a 
tumour-suppressor gene, DLC1, 23 it was of interest to determine 
its expression in genome-modified CLECs. We performed 
quantitative reverse transcription (RT)-PCR to compare DLC1 
transcript levels (exons 1-2 and exons 8-9) in nine 8p22- 
integrated oligoclones and in unmodified control CLECs from 
the same donor (that is, electroporated without plasmid DNA but 
having the same number of population doublings in culture). 
DLC1 transcript levels were normalized to glyceraldehyde-3- 
phosphate dehydrogenase (GAPDH) expression in all reactions. 
The ratio of DLC1 transcript levels in CLECs with transgene 
integration in 8p22 relative to DLC1 transcripts levels of control 
CLECs ranged between 0.58-1.603 (exons 1-2) and 0.851-1.38 
(exons 8-9), showing minimal perturbations to DLC1 expression 
caused by transgene integration at 8p22 (Figure 3a). 

Transcriptome and genome copy number analyses of 
DLC1 -integrated oligoclonal CLECs 

The extent to which transgene integration in DLC1 altered global 
gene expression of oligoclonal CLECs was investigated by 
transcriptome profiling on Affymetrix PrimeView expression arrays 
(Affymetrix, Santa Clara, CA, USA). Transcripts that differed by ^ 2.5- 
fold compared with unmodified control CLECs were identified. 



Transcriptome analysis identified a total of 341 genes out of a 
total of 17 719 expressed genes (1.9%) that were commonly altered 
in expression in all 8 oligoclonal CLECs with DLC1 integration 
compared with control CLECs (Supplementary Table 2). Of these, 93 
genes (0.5%) were upregulated and 248 genes (1.4%) were 
downregulated. No changes were detected in DLC1 expression or 
genes that mapped within a 1-Mb window (LONRF1, KIAA456 and 
C8orf48) centred on the 8p22 integration site for all eight oligoclonal 
CLECs. Genes whose expressions were altered in oligoclonal CLECs 
were classified and clustered according to their biological functions 
and their proportional categorizations are depicted in Figures 3b 
and c. It was noteworthy that up to 40% of underexpressed genes 
were involved in cell cycle regulation, suggesting that oligoclonal 
CLECs with DLC1 integration could be less proliferative compared 
with control CLECs. A possible reason for downregulated cell cycle 
genes was the increased metabolic burden from artificially induced 
overexpression of a transgenic protein. 24 Given that oligoclonal 
CLECs had been cultured continuously in vitro for at least 2 months, 
downregulation of cell cycle genes could also reflect onset of 
senescence of these primary cells. In either case, downregulation of 
cell cycle genes was contrary to the transcriptional profile of 
transformed cells. Further characterization of oligoclonal CLECs by 
the MTS assay revealed decreases in proliferation in five of six 
oligoclonal CLECs evaluated (15, 39, 48, 50, 59 and 62) by 
28% -47% (Supplementary Figure S3). This was in accord with 
the observed downregulated transcription of cell cycle genes that 
caused a slowdown in proliferation of these cells. 

The list of genes with altered expression was submitted to 
DAVID pathway mapping analysis and compared against a list of 
known and potential oncogenes. 25 This clustered downregulated 
genes (46 genes) mainly to cell cycle and DNA repair pathways 
(KEGG pathway database; Supplementary Figure S2). Genes with 
upregulated expression were not significantly clustered to any 
pathway. Referencing the altered genes to a database of known 
and potential oncogenes identified 9 and 31 oncogenes that were 
up- and downregulated, respectively (Supplementary Table 2). 
In summary, transcriptome analysis of oligoclonal CLECs with 
transgene integration in intron 7 of DLC1 did not display 
signatures of transformed cells. 

With a single exception, all eight oligoclones with DLC1 
integrations did not show significant alterations in genome copy 
number determined by hybridization on high-density single- 
nucleotide polymorphism arrays (average probe spacing 
11 48 bp; Affymetrix CytoScan HD array). A single amplification 
(copy number state of 3) in the peri-centromeric region of 
chromosome 1 9 (spanning 41 1 8 kb) was present in oligoclone 1 0 
(Figure 4). Only a single gene, ZNF254, resided within this 
amplified region but its expression was not changed by more 
than two-fold. 



► 

Figure 2. Identification of oligoclonal CLECs with transgene integration at 8p22 and FVIII secretion, (a) Oligoclonal cells (four flow-sorted cells 
per well) from a bulk population of G418-selected CLECs electroporated with paffBhybrid FVIII and phiC31 integrase were investigated by 
direct in situ PCR for transgene integration at the 8p22 hot spot. In situ lysed cells were screened using Phusion human specimen direct PCR kit 
(Thermo Scientific) and primers specific for vector and genomic sequences at the 8p22 integration site to detect the presence of left and right 
integration junctions. Intron-spanning vector-specific primers were used to amplify C1 domain of the integrated FVIII cDNA (FVIII vector). 
Control genomic PCR amplified a 900-bp region in 19q 13.42 {AAVS1 locus), -ve' denotes minus template PCR amplification while ' + ve' 
denotes amplification from genomic DNA isolated from a previously identified clonal CLEC with integration at 8p22. Amplified products were 
electrophoresed on 1% agarose gels and imaged using BioRadGel Doc 2000 transilluminator and QuantityOne software. Red box highlights 
samples, which were positive for left and right integration junctions, FVIII vector and control PCR. (b) Genome-modified oligoclonal CLECs 
identified by junction PCR to be positive for transgene integration at chromosome 8p22 were screened by FISH. Fluorescence images of 4,6- 
diamino-2-phenylindole (DAPI)-stained cells hybridized with fluorescein isothiocyanate (FITC)-labelled probes specific to integrated vector 
(green signal) and Texas Red-labelled centromeric probes specific to chromosome 8 (red signal) are shown (original magnification x 600). 
Integration of transgene at chromosome 8 is indicated by the close proximity of the vector-specific green signal and chromosome 8 
centromere red signal, (c) CLECs that were electroporated without any plasmid DNA (EP only) and oligoclonal CLECs that were positive for 
transgene integration at 8p22 (10, 15, 16, 28, 39, 40, 50, 59, 62, 68 and 69) and three oligoclonal populations that were negative for 8p22 
transgene integration (18, 30 and 47) were evaluated for FVIII secretion (day 40 post-electroporation). Overnight conditioned media was 
assayed for FVIII activity using a Coamatic FVIII kit (Chromogenix). Data are mean±s.e.m.; n = 3 per group. 
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Figure 3. RT-PCR and transcriptome analysis of 8p22 oligoclonal 
CLECs, (a) Quantitative RT-PCR was performed on control CLECs that 
had been electroporated without any DNA and oligoclonal CLECs 
with transgene integration at 8p22 to quantify DLCI transcript levels 
(exons 1-2 and exons 8-9) and GAPDH expression. The mean fold 
change in G/APDH-normalized DLCI expression levels in oligoclonal 
samples (exons 1-2 and exons 8-9) are shown relative to that of 
control CLECs as determined by the 2" AACt method. Data are 
mean ± s.e.m.; n = 3 experiments per group, (b, c) Categorization of 
genes that differed in their expression by ^ 2.5-fold in all eight 
oligoclonal CLECs relative to control CLECs. The percentage of (b) 
significantly upregulated and (c) significantly downregulated genes 
categorized according to their biological functions are depicted. 



Tumourigenicity and FVIII expression in vivo 
FVIII-expressing DZ_C7-integrated oligoclonal CLECs were 
implanted in the subcutaneous nuchal region of NOD-SCID IL2Ry 
null (NSG) mice. Survival and engraftment of implanted cells was 
monitored by measuring levels of transgenic hybrid FVIII secreted 
by implanted CLECs into murine plasma using a FVIII-antigen 
capture assay specific for human FVIII that did not cross-react with 
murine FVIII. Untreated mice and mice implanted with CLECs that 
were electroporated without plasmid DNA served as negative 
controls for tumour formation and plasma hybrid FVIII measure- 
ments, respectively. NSG mice implanted with Hs746T, a human 
cancer cell line, served as positive tumourigenic controls. 



Hybrid FVIII was readily detected at high levels in murine 
plasma 3 days post-implantation (335.65 ± 19.98 mU ml -1 ; 
Figure 5a). Although plasma levels of transgenic FVIII showed a 
declining trend (that is, 236.9 ± 24.67 mU ml -1 on day 7 and 
40.1 7 ± 6.29 mU ml -1 on day 14), day 14 levels remained 
significantly higher than untreated mice (10.5 ± 2.45 mU ml" 1 ; 
P = 0.0094) and mice implanted with unmodified CLECs 
(8.67 ± 1.06 mU ml" 1 ; P = 0.0063). Plasma FVIII levels were not 
different from control mice by day 30 (P>0.05) although CLECs 
were demonstrable at implantation sites for at least 3 months 
after implantation (Figure 5b). Loss of transgenic FVIII levels in 
plasma may not have been caused by silencing of FVIII transgene 
expression in implanted CLECs. We recovered and re-established 
secondary cultures from implanted cells, and documented that 
they retained the capacity to secrete FVIII in vitro (1 304 ± 34.05 mU 
FVIII/10 6 cells per 24 h). Recovered CLECs were positive for human 
vimentin by immunohistochemical staining (Figure 5c). Collec- 
tively, these data indicate that a proportion of the initially 
implanted cells did survive and established long-term engraft- 
ment but were not sufficiently numerous to secrete detectable 
levels of transgenic FVIII in murine plasma. Although no tumours 
were observed for up to 3 months after implantation of 3 x 10 6 
unmodified {n = 4) or FVIII-expressing CLECs (n = 16), large 
tumours were readily observed as early as 3 weeks in mice 
implanted with the same number of Hs746T cells (n = 4; 
Figure 5d). These data showed that oligoclonal DLC 7-integrated 
CLECs, which had been cultured in vitro for at least 2 months were 
unlikely to be tumourigenic. 

DLC1 integration in primary adult cells 

Primary cultures of bone marrow-derived stromal cells, adipose- 
derived stromal cells and dermal fibroblasts were co-electropo- 
rated with pctffShybrid FVIII and phiC31 integrase plasmids, then 
selected with G418 to derive stable cell cultures. Integration 
junction PCR was positive for site-specific transgene integration in 
intron 7 of DLC1 in all three primary cell types (Figure 6a), showing 
that this site was permissive for transgene integration in a range 
of adult primary human cell types. Moreover, these stably 
modified cells secreted FVIII expression for at least 1 month 
(Figure 6b) indicating that the human ferritin light chain promoter 
was not silenced and could stably drive hybrid FVIII expression in 
several cell types. 



DISCUSSION 

Oncogenic complications that have marred several clinical cell 
therapy trials have renewed concerns about the risks of randomly 
integrating gene delivery vectors. Developing non-viral methods 
capable of modifying genomes efficiently and with precision will 
be necessary to advance the biosafety of gene-based cell 
therapies. 

Although phiC31 integrase is effective in mediating stable 
transgene integrations and expression in a variety of mammalian 
cells, its activity in integrating donor DNA into multiple genomic 
sites could be perceived to be a barrier to clinical applications. 
Along with others, 12 ' 13 ' 20 we 17 have noted that intron 7 of DLC 1 in 
8p22 is a frequent site for phiC31 integrase-mediated transgene 
insertion. Its permissiveness to genome modifications is possibly 
due to wide expression of DLC 1 in normal human tissues. 26,27 Here 
we show that the use of phiC31 integrase for stable transgene 
integration could still be useful and rendered safer by screening, 
identifying and evaluating clonal cells that have integrated the 
transgene into a specific site, that is, DLCi. Our data suggest that, 
given the high frequency of integration in DLCi, such cells could 
justifiably be evaluated for insertional oncogenesis and efficacy of 
transgene expression ex vivo by multiple techniques to determine 
their suitability for cell therapy. DLCi is not an obvious safe 
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genomic harbour site, being a tumour suppressor. Moreover, it 
was not known if this site is capable of sustaining durable 
transgene expression. It was therefore of interest to assess the 
genotype and phenotype of primary human cells that had 
integrated a hybrid FVIII transgene within DLC1. We have 
previously shown that transgene integrations into 8p22 were 
either mono- or biallelic. 17 Hence, the foregoing genomic and 
transcriptome analyses in this study reflected genome 
modification at this single site, unconfounded by the effects of 
simultaneous transgene integrations into heterogeneous sites. We 
used integration junction PCR as a screening method to identify 
cells with transgene integration at 8p22, followed by FISH with 
transgene-specific probes to confirm integration in chromosome 8. 
Although these techniques were facile as screening methods for 
single site-specific integration, Southern blotting could have 
detected the occurrence of transgene integrations at multiple 
genomic sites. 

Some oligoclones (#55, 20, 21, 24 and 30) appeared to be 
positive only for the left integration junction and not the right 
integration junction. In principle, deletion of vector and/or 
genomic DNA sequences where the integration junction primers 
were designed to anneal would cause the integration junction 
PCR to fail. Deletions of vector sequences as well as genomic DNA 
following phiC31 integrase-mediated integration have been noted 
in our current study (Supplementary Figure Sic) and in previous 
studies. 17,29 The vector-specific and genome-specific primer 
sequences for the right integration junction were 37 and 264 bp, 
respectively, from the integration (cross-over) region. Deletions to 
the vector of >37bp or to the genomic DNA of > 264 bp would 
thus have deleted the primer binding sites, resulting in negative 
amplification of the right integration junction. 

We showed by RT-PCR and transcriptome analyses that DLO 
expression in nine oligoclonal transgene-integrateo! CLEC cultures 



was 0.6-1.6-fold relative to naive control CLECs. Insertional 
oncogenesis events that complicated recent clinical studies have 
been ascribed to transactivation of neighbouring oncogenes such 
as LM02, CCND2, BMI1 and the MDS-EVI1 complex. 7 " 9 Our 
transcriptome data show that intronic integration had little 
effect on DLO expression but also did not dysregulate 
transcription of neighbouring genes within 1 Mb centred on the 
integration site. Moreover, transcriptome profiles of oligoclonal 
CLECs with intragenic DLO transgene integration did not bear the 
hallmarks of transformed cell lines as a large majority of cell cycle 
genes were downregulated, suggesting that these genome- 
modified cells were likely to be less proliferative than 
transformed cells. Genome-wide copy number profiling of 
oligoclonal CLECs on a high-resolution single-nucleotide 
polymorphism array revealed no significant alteration in seven 
of eight oligoclonal CLECs. Only a minor amplification of a 41 18-kb 
peri-centromeric region in chromosome 19 was observed in a 
single oligoclone that did not alter transcription of the only gene 
in close proximity, ZNF254. Given the high-density coverage of 
probes used in the array, the absence of significant copy number 
changes suggests that intragenic DLO integration did not induce 
structural rearrangements of chromosomes such as translocations 
as these are usually unbalanced in nature. Although structurally 
abnormal chromosomes have been reported with other phiC31 
integrase studies, 28,29 the karyotypes reported were not selective 
for cells with 8p22 transgene integration. 

The benign safety profile of phiC31 integrase-modified oligo- 
clonal CLECs suggested by in vitro analyses was confirmed by 
absence of in vivo tumourigenicity in NSG mice. Although no 
tumours were observed for at least 4 months post-implantation of 
3 million transgene-modified oligoclonal CLECs, parallel implanta- 
tion of 3 million of a tumourigenic carcinoma cell line (Hs746T) 
resulted in tumour formation within 3 weeks. Given that tumours 
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Figure 5. Assessment of in vivo FVlll secretion and tumourigenicity following implantation of CLECs stably secreting transgenic hybrid FVlll in 
NSG mice, (a) Plasma FVlll levels of NSG mice that were unimplanted (n = 4) or implanted subcutaneously with 3 x 10 6 cells per animal of 
(i) Matrigel-encapsulated CLECs that received electroporation only (EP only, n = 4); or (ii) FVIII-secreting oligoclonal CLECs (n = 16) with 
paffBhybrid FVlll integrated at 8p22 were measured using an assay specific for human FVlll (Visulize FVlll antigen ELISA kit) on days 3, 7, 14 
and 30 post-implantation. Plasma FVlll levels of mice implanted with FVIII-secreting cells were significantly higher (P<0.05, Student's unpaired 
f-test) on days 3 (yellow bar), 7 (orange bar) and 14 (grey bar) compared with control mice that were either unimplanted or implanted with 
cells electroporated without any plasmid DNA (EP only). Data are mean ± s.e.m. (b) Immunohistochemical staining of Matrigel implants 
retrieved by excision from implanted mice at the indicated time points post-implantation show evidence (brown staining) for the presence of 
engrafted CLECs expressing human vimentin. Scale bar = 50 jam. (c) (Left) Brightfield view of outgrowth cells from a culture of Matrigel 
implants retrieved 2 months post-implantation; and (right) immunohistochemical staining of explant cultured cells for human vimentin 
expression. Scale bar = 1 00 jam. (d) No tumours were observed in mice implanted with Matrigel-encapsulated 3 x 1 0 6 genome-modified CLECs 
at day 90 post-implantation (left) while nodular tumours developed in mice implanted with equal numbers of Matrigel-encapsulated Hs746T, 
a known tumourigenic cell line, 1 -month post-implantation (right). Arrow points to Matrigel implant. 



formed in this strain of mice implanted with as few as 100-1000 
cancer cells, 30,31 even small numbers of genome-modified CLECs 
would have generated tumours if they had been oncogenically 
transformed by transgene integration in DLC1. Taken together, our 
biosafety assessments indicated that intragenic DLC1 transgene 
integration did not predict insertional oncogenicity, and genome- 
modified CLECs were non-tumourigenic in vivo. The eventual 
decline of plasma FVlll levels to basal levels by day 30 may have 
resulted from loss of a substantial proportion of implanted cells 
because of sub-optimal engraftment and vascularization in a 
xenogeneic host. Transgene silencing in vivo was another possible 
cause. We noted that although secondary cultures established 
from explanted tissues still secreted FVlll in vitro (1304 ±34.05 mil 
FVIII/1 0 6 cells per 24 h) following engraftment in mice for at least 3 
months, these levels were lower than pre-implantation levels 
(261 0-5724 mU FVIII/1 0 6 cells per 24 h), suggesting some level of 
transgene silencing following implantation. Collectively, these 
data indicate that a proportion of the initially implanted cells did 
survive and established long-term engraftment but were not 



sufficiently numerous to secrete detectable levels of transgenic 
FVlll in murine plasma. A comprehensive evaluation of durable 
FVlll expression in vivo may require long-term engraftment of 
genome-modified cells in an autologous model. 

Unlike semi-random transgene integrations mediated by viral 
vectors, several non-viral methods using meganucleases, zinc- 
finger nucleases, 33 nickases, 34 transcription activator-like effector 
nucleases 35 and the CRISPR/Cas9 system 36 have been developed 
recently to enable site-directed integration into a pre-determined 
genomic locus. The choice of genomic sites that are predictably 
safe from insertional oncogenesis and that will not silence the 
integrated transgene remains challenging. Thus far, AAVS1 
(19q13.42) and CCR5 (3p21 .31) are the two most common 
genomic loci selected for site-directed integrations. 10 Of these, 
integration into the AAVS1 locus is thought to be benign although 
there is as yet no comprehensive determination of its safety. This 
locus, which falls within the first intron of PPP1R12C, supports 
robust transgene expression in several cell types. However, as the 
function of the PPP1R12C protein is largely unknown, cellular and 
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Figure 6. Transgene integration at 8p22 and durable FVIII secretion 
in phiC31 integrase-modified primary adult human cells, (a) G418- 
resistant stable cells derived from primary cultures of bone marrow- 
derived stromal cells (BMSC), adipose-derived stromal cells (ADSC) 
and normal dermal fibroblasts (NF123) co-electroporated with 
paffBhybrid FVIII and phiC31 integrase were examined for evidence 
of transgene integration into intron 7 of DLC1 by left and right 
integration junction PCR. -ve' refers to minus template amplifica- 
tion. Left and right integration junctions were amplified with vector- 
specific and genomic DNA-specific primers while control genomic 
PCR was performed with a pair of genome specific-primers 
amplifying a 900-bp region within chromosome 19. Amplified 
products were electrophoresed on 1% agarose gels and imaged 
using BioRad Gel Doc 2000 transilluminator and QuantityOne 
software. A composite image of the original gel is shown. Black 
arrows indicate the predicted integration junction PCR amplified 
bands, (b) FVIII activity in overnight conditioned media of untreated 
(Wt) and G418-resistant (phiC31 integrase-modified) BMSC, ADSC 
and NF123, determined 1 month after co-electroporation with 
paffBhybrid FVIII and phiC31 integrase. Data are mean ± s.e.m.; n = 3 
per group. 



organismal consequences of its altered expression remain to be 
ascertained. Integration in AAVS1, CCR5 and DLC1 sites support 
transgene expression probably because they are intragenic and 
transcriptionally active in many cell types, unlike intergenic sites 
that may be heterochromatic or otherwise unpermissive for 
transgene expression. 37 Intragenic transgenes are also less likely 
to perturb intergenic regions whose functions have yet to be 
identified and understood. Validation of candidate safe genomic 
harbours is urgently needed for gene and cell therapy to progress 
to clinical practice with acceptable risk-benefit calculus. Our 
characterization of transgene integration into DLC1 has several 
features that could be expected of a genomic safe harbour. 
Moreover, we show that DLC1 integration supports durable 
transgene expression comparable to other putative genomic safe 
harbours {AAVS1 and CCR5). It resulted in minimal changes to 
the transcriptome and genome copy number, and CLECs with 
DLC 7-integrated transgene were not tumourigenic in vivo. This 
was not surprising because although DLC1 is a tumour suppressor, 
its expression was not extinguished by the integrated transgene. 
However, in the event that DLC1 expression could decrease 
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significantly in other situations, concurrent expression of 
the closely related DLC2 (13q12), a RhoGAP with similar specificity 
for RhoA and Cdc42, 38 may functionally compensate for reduced 
activity of DLCh 

Increasing the number of genomic sites into which transgenes 
can be integrated without insertional oncogenesis will be helpful 
for the potentially wide range of gene and cell therapy 
applications. The same site may well be selectively oncogenic 
in some, but not all, cell types because gene products often 
have cell type-specific functions. Intronic integration that does 
not change the expression of the host gene, as in our DLC1 data, 
is doubly advantageous as it is more likely to be functionally 
benign yet permissive of transgene integration and durable 
expression. 

Our approach suggests that integrating a therapeutic transgene 
in primary human cells ex vivo with precision at a preselected 
single genomic site should mitigate risks of direct in vivo gene 
delivery by enabling comprehensive characterization of the 
transcriptome, genome copy number, karyotype and cellular 
behaviour. This strategy is feasible in primary cell types that are 
simple to procure, expandable in vitro to clinically meaningful 
numbers and genomically stable. Clonal cultures of genome- 
modified cells could also be more readily derived from primary 
cells having these characteristics. However, oligoclonal CLECs 
cultured for about 2 months had downregulated transcription of 
cell cycle genes and reduced proliferation. This development can 
be anticipated also for other primary human cells such as 
fibroblasts, tissue-specific progenitor and stromal cells after 
long-term culture. Nonetheless, finite expandability of CLECs does 
not necessarily preclude their application for haemophilia therapy. 
Estimates based on FVIII expression from transgene-modified 
CLECs (2131 mU FVIII/10 6 cells per 24 h) suggest that approxi- 
mately 100-400 million cells are required to achieve a steady-state 
plasma FVIII level of 10% in a 20-40 kg child. Given that 
approximately 20 million transgene-modified CLECs can be 
derived from each oligoclonal CLEC identified by our screening 
procedure, this can be achieved by identifying 5-20 oligoclonal 
CLECs to generate the cell numbers required for clinically effective 
cell therapy. Based on our current study, identification of 5-20 
oliogoclonal CLECs by screening for transgene integration at 8p22 
is technically feasible. These clonal cells can be readily expanded 
to 100-400 million cells in total before growth begins to slow 
down. Alternatively, induced pluripotent stem cells that have 
unlimited proliferation capacity may be better suited in practice 
for transgene integration and clonal selection of cells with site- 
specific integration at 8p22, provided concerns pertaining to 
tumourigenicity are addressed. 

We chose haemophilia A as a model disease caused by 
deficiency of a single protein to demonstrate the feasibility of 
autologous cell therapy using genome-modified CLECs. Secretion 
of BDD-hybrid FVIII was at least fivefold higher than BDD-human 
FVIII. Given the inherently poor expression of recombinant human 
FVIII, 39 " 41 the more efficiently secreted BDD-hybrid FVIII 42 could be 
the preferred transgene in future preclinical studies as fewer cells 
would be required to achieve therapeutically significant plasma 
FVIII levels. 

The convergence of new methods for precision genome 
modification, high-throughput techniques for characterizing 
whole genomes and transcriptomes and improved design of 
donor transgenes could advance gene-based cell therapy of 
inherited and acquired disorders among paediatric and adult 
patients. 



MATERIALS AND METHODS 

Hybrid FVIII plasmid assembly 

The hybrid porcine-human FVIII construct consisted of porcine Al and A3 
domains, human signal peptide, A2, and partial B domains (retaining the 
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C-terminal 266 amino acids with eight glycosylation sites), C1 and C2 
domains. Human FVIII components were derived by PCR amplification from 
the full-length sequence in pSP64-F8 (American Type Culture Collection, 
Manassas, VA, USA). Porcine FVIII domains were derived by RT-PCR of 
porcine liver RNA. An overlap PCR strategy was used to assemble 
components of the hybrid FVIII cDNA. The construct, expressed from the 
human ferritin light chain promoter, was eventually cloned in a vector 
containing the 350-bp attB fragment 22 and neomycin resistance gene. The 
resulting construct was named patt£hybrid FVIII. 



Cell culture 

Primary human dermal fibroblasts (KF1 and NF123), human adipose- 
derived stromal cells, human bone marrow-derived stromal cells and CLECs 
were provided by CellResearch Corporation, Singapore. Dermal fibroblasts, 
adipose-derived stromal cells and bone marrow-derived stromal cells were 
cultured in Dulbecco's modified Eagle's medium and 25 itim glucose 
supplemented with 10% foetal bovine serum. The isolation and culture of 
CLECs has been described previously. 1 7,43,44 

Gene transfer and stable cell generation 

All gene transfer studies were performed using a 4D-Nucleofector device 
(Lonza). CLECs were electroporated with Nucleofector solution P1 and 
setting CM113, bone marrow-derived stromal cells and adipose-derived 
stromal cells with Nucleofector solution P2 and setting EW104 while 
dermal fibroblasts (NF123) with Nucleofector solution P2 and setting 
FF113. Transfection efficiency was gauged by fluorescence-activated cell 
sorting analysis for green fluorescent protein-positive cells 24 h post- 
electroporation with 5 |ig of pmaxGFP vector (Lonza). Cells stably modified 
by phiC31-integrase were generated by co-electroporation with 12|ig of 
pattBhybrid FVIII and 2.5 |ig pCMV-Int plasmid 22 followed by selection with 
1 mg ml -1 of G418 for 5 days starting from day 6 post-electroporation. 

Integration junction PCR and sequencing 

Integration junction PCR was performed on 200 ng of genomic DNA (DNA 
extraction kit; Favorgen Biotech Corp., Ping-Tung, Taiwan) using DyNAzyme 
EXT DNA polymerase (Thermo Scientific, Rockford, IL, USA). Left integration 
junction PCR (454-bp amplicon) was performed using primer pairs specific to 
a locus on chromosome 8p22 (forward: S'-GGGCTCTGGAGTAAAGGTGAAA-^) 
and donor vector (reverse: S'-GTTCGCCGGGATCAACTACC-^). Right integra- 
tion junction PCR (333-bp amplicon) was performed using primer pairs 
specific to the vector sequence (forward: 5'-TCGACGATGTAGGTCACGG-3') 
and chromosome 8p22 genomic DNA (reverse: S'-GCATGGCCTCATTTC 
CGTCT-3'). Transgene-specific PCR was performed with intron-spanning 
primers specific to the CI -domain of the integrated hybrid FVIII cDNA 
(forward: 5'-AGCnTCAGAAGAGMCCCGACAC-3'; reverse: 5'-TCCCAGGGGAG 
TCTG AC ACTTCTTG CTGTAC ACC AGG A AAGT3' ) . Control genomic PCR (900-bp 
amplicon) was performed with genomic primers (forward: 5'-AAGAAGCG 
CACCACCTCCAGGTTCT-3'; reverse: 5 , -ATGACCrCATGCTCTTGGCCCrCGTA-3 / ). 
Amplified products were electrophoresed on 1% agarose gels and imaged 
using BioRad Gel Doc 2000 transilluminator (Bio-Rad Laboratories, 
Hercules, CA, USA). 

Screening of oligoclonal CLECs with transgene integration at 8p22 
To screen oligoclonal CLECs for transgene integration at chromosome 
8p22 locus, 45 G418-resistant stable population of CLECs (pafffihybrid 
FVIII + pCMV-Int, G418 selected) were flow sorted (four cells per well) 
into individual wells of 96-well plates and allowed to expand in culture. 
A replica plate for further analysis of clones was established by detaching 
and re-seeding cells into a second 96-well plate. Cells for 8p22 integration 
screening (in 96-well plates) were lysed in situ with 60 |il of lysis buffer and 
screened by direct PCR for the presence of integration junctions and the 
integrated transgene using Phusion Human Specimen Direct PCR kit 
(Thermo Scientific) and the same primers above. 

Fluorescence in situ hybridization 

FISH to detect integration of partflhybrid FVIII was performed on 
interphase nuclei of genome-modified CLECs with fluorescein-12-dUTP 
(PerkinElmer, Waltham, MA, USA) random prime labelled probe generated 
by PCR amplification of neomycin cDNA from pattB hybrid FVIII 
(817-bp PCR product: forward primer S'-TTGCACGCAGGTTCTCCGGC-S'; 
reverse primer S'-GGCGTCGCTTGGTCGGTCAT-S'). Texas-Red-5-dUTP-labelled 



human chromosome 8 centromeric probe (Children's Hospital Oakland 
Research Institute, Oakland, CA, USA) was used to determine colocalization 
with chromosome 8. Hybridized slides were counterstained with 4,6- 
diamino-2-phenylindole and visually enumerated for probe signals with an 
Olympus BX61 epifluorescence microscope (Olympus, Tokyo, Japan). 



FVIII assay 

Factor VIII activity in conditioned media was measured using a 
Chromogenix Coamatic Factor VIII kit whose detection limit was 1% FVIII 
activity (Chromogenix, Molndal, Sweden). Human-specific FVIII antigen 
levels in murine plasma was quantified using Visulize FVIII antigen ELISA kit 
having a detection limit of 0.8% Factor VIII antigen (Affinity Biologicals, 
Ancaster, ON, Canada). Both assays were performed according to the 
manufacturers' protocols. 



Transcriptome analysis 

Total RNA isolated from control CLECs (electroporated without plasmid 
DNA and having the same number of population doublings) and genome- 
modified oligoclonal CLECs with integrations at 8p22 served as starting 
material for transcript profiling on GeneChip PrimeView Human Gene 
Expression array (Affymetrix) following the recommended protocol. 
Transcription expression data were analyzed using GeneChip Operating 
Software (Affymetrix). Transcripts whose expression levels differed by > 2.5- 
fold in genome-modified oligoclonal CLECs compared with naive CLECs were 
considered significantly altered and were further analyzed. DAVID (Database 
for Annotation, Visualization and Integrated Discovery) 2.1 Functional 
Annotation Tool (http://david.abcc.ncifcrf.gov) was used to ascribe functions 
and other annotations for significantly altered transcripts and for pathway 
mapping. Altered transcripts were also referenced to a compilation of known 
proto-oncogenes and tumour-suppressor genes (http://www.bushmanlab. 
org/links/genelists). 



Genome copy number change analysis 

Genome-modified oligoclonal CLECs with 8p22 integrations and control 
CLECs (electroporated without plasmid DNA and having the same number 
of population doublings) were evaluated using Cytoscan HD array and 
Chromosome Analysis Suite for data analysis (Affymetrix). Regions of copy 
number gain or loss were defined as segments of ^50 consecutive probes 
concordant for significant copy number abnormalities. Log 2 signal 
intensity ratios >0.3 and ^-0.3 were criteria for significant copy number 
gain and loss, respectively. 



RT-PCR analysis 

RT-quantitative PCR was used to determine changes in DLC1 transcript 
levels in clonal CLECs with transgene integration in 8p22. CLECs 
electroporated without plasmid DNA and having the same number of 
population doublings were controls. In brief, total RNA was extracted using 
TriPure isolation reagent (Roche Applied Science, Indianapolis, IN, USA), 
treated with DNase I (Qiagen, Hilden, Germany) to remove contaminating 
genomic DNA and purified using RNeasy Mini kit (Qiagen). Reverse 
transcription was performed on approximately 1 |ig of DNAse-treated RNA 
using iScript Advanced cDNA synthesis kit (Bio-Rad Laboratories), 
according to the manufacturer's instructions. Quantitative-PCR was 
performed using 2|il of first strand cDNA in a 20|il reaction volume using 
GoTaq qPCR Master Mix (Promega Corp., Madison, Wl, USA) and 45 cycles 
at an annealing and extension temperature of 62 °C using CFX96 Real-Time 
PCR detection system (Bio-Rad Laboratories). Intron-spanning exonic 
primers were used to amplify DLC1 exons 1-2 (forward: 5'-TCCTG 
CCCCAATGGAATGTC-3'; reverse: 5'-GTTGGTGTGCCTGATGGAGA-3'), exons 
8-9 (forward: 5'-GAAGGGGATGCAGCGGATAG-3'; reverse: 5'-AGCAGGG 
CCGTTAGCTTTAG-3') and a housekeeping gene, GAPDH exon 6-7 (forward: 
S'-GCCTCCTGCACCACCAACT-S'; reverse: S'-CGCCTGCTTCACCACCTTC-S'). 
DLC1 transcript levels were normalized to GAPDH expression levels and 
the fold-change in DLC1 transcript levels in clonal CLECs with 8p22 
integration was reported relative to DLC1 transcript levels in control CLECs, 
using the 'delta-delta Ct method'. 46 

Cell proliferation assay 

Oligoclonal CLECs that had been in culture for at least 2 months and 
control CLECs that were electroporated only and cultured for only 1 month 
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were seeded at 1000 cells per well and allowed to proliferate for 3 days, 
after which MTS assay (CellTiter 96 Aqueous One Solution Cell Proliferation 
Assay, Promega Corp.) was performed to quantify live cells by absorbance 
readings at 490 nm. 

Mice and implantation studies 

All animal studies were conducted in accordance with ethical and 
international animal use guidelines following approval by SingHealth 
Institutional Animal Care and Use Committee (IACUC). The tumourigenic 
potential of hybrid FVIII-expressing clonal CLECs with transgene integra- 
tion at 8p22 locus was investigated by implanting 3x10 6 Matrigel- 
encapsulated (BD Biosciences, San Jose, CA, USA) transgenic cells into the 
subcutaneous nuchal region of NOD-SCID IL2Ry null (NSG) mice (Jackson 
Laboratories, Bar Harbor, ME, USA; n = 16). Untreated mice (n = 4) and 
mice implanted with CLECs that received electroporation only (n = 4) 
served as negative controls for tumour formation and plasma hybrid FVIII 
measurements, respectively. NSG mice implanted with a tumourigenic cell 
line (Hs746T, American Type Tissue Culture, Manassas, VA, USA) served as 
positive tumourigenic controls. Survival and engraftment of implanted 
cells was monitored by measuring levels of transgenic hybrid FVIII 
(secreted by implanted CLECs) in murine plasma using a FVIII-antigen 
capture ELISA assay (Affinity Biologicals). Mice were monitored weekly for 
visible tumours at implantation sites. Matrigel implants were excised under 
sterile conditions to re-establish secondary cultures and for immunohis- 
tochemistry. Secondary explants cultures were established by dicing 
excised Matrigel implants into 0.5 mm cubes, seeding them onto 10-cm 
tissue culture dishes and culturing them exactly as for primary CLECs. 

Immunohistochemistry 

Immunohistochemical staining to identify implanted human cells was 
performed on formalin-fixed paraffin-embedded tissues excised from 
implantation sites. Five-micron tissue sections blocked in 10% goat serum 
for 30min were incubated with 1:3 diluted mouse anti-human vimentin 
antibody (Clone V9, Zymed, San Francisco, CA, USA) for 1h at 37 °C 
followed by 30-min incubation with 1:2 diluted DAKO REAL EnVision 
rabbit/mouse horseradish peroxidase-conjugated secondary antibody 
(Dako, Glostrup, Denmark). Tissue sections were washed three times with 
phosphate- buffered saline between each incubation step. Slides were 
stained with DAB chromogen (1 part DAB chromogen + 3 parts DAKO 
REAL substrate buffer) for 3 min, counterstained with Mayer's haematoxylin 
and mounted with DPX mount solution. Tissue sections were visualized 
and imaged with an inverted microscope (Axiovert 25CFC, Carl Zeiss, 
Oberkochen, Germany) under x 20 magnification using the KS400 
software (Carl Zeiss). 

Statistical analyses 

Results are expressed as mean ± s.e.m. or s.d. where appropriate, with 'n' 
indicating the total number of samples or measurements. Student's 
unpaired f-test with two-tailed P-values and 95% confidence interval was 
used for comparison between two groups. P<0.05 was considered 
statistically significant. 
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